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The generation and control of quantum states of light constitute fundamental tasks in cav- 
ity quantum electrodynamics (QED)EM^. The superconducting realization of cavity QED, 
circuit 

qEdEQElI enables on-chip microwave photonics, where superconducting qubits ^^"^ 
control and measure individual photon states^HSD. a long-standing issue in cavity QED 
is the coherent transfer of photons between two or more resonators. Here, we use circuit 
QED to implement a three-resonator architecture on a single chip, where the resonators 
are interconnected by two superconducting phase qubits. We use this circuit to shuffle 
one- and two-photon Fock states between the three resonators, and demonstrate qubit- 
mediated vacuum Rabi swaps between two resonators. This illustrates the potential for 
using multi-resonator circuits as photon quantum registries and for creating multipartite 
entanglement between delocalized bosonic modest 

The combination of high-finesse electromagnetic cavities with atoms or qubits enables fundamental 
studies of the interaction between light and matter. The cavity provides a protected environment for stor- 
ing and tailoring individual photonic excitations'^^. Both stationary^!! an d propagating^^ nonclas- 
sical fields can be synthesized using such systems, enabling quantum memory and quantum messaging 8 . 
A critical challenge however is the extension from single to more versatile multi-cavity architectures* 27 * 29 -^, 
allowing manipulation of spatially separated bosonic modes. While the entanglement of different modes 
of a single cavity has been shown in atomic systems^!! and a coupled low- and high-quality factor res- 
onator studied in circuit QED^ 2 ^, coherent dynamics between two or more high-quality factor cavities 
has yet to be demonstrated. Here we describe a triple-resonator system, where three high-quality factor 
microwave resonators are coupled to two superconducting phase qubits (cf. Fig. [T]). The qubits serve as 
quantum transducers that create and transfer photonic states between the resonators^. The quantum 
transduction is performed by means of purely resonant qubit-resonator interactions, rather than disper- 
sive coupling 27 , enabling rapid transfers between resonators with significantly different frequencies. As 
an important example, we demonstrate single-photon Rabi swaps between two resonators de-tuned by 
~ 12000 resonator linewidths. 

Figure []Ji,b shows the main experimental elements, which comprise three coplanar waveguide res- 
onators, R a , and R c , two phase qubits, Q x and Q 2 , and two superconducting quantum interference 
devices (SQUIDs) used for qubit state readout. Each qubit is coupled to a control line which is used to 
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Fi gure ll Experimental architecture and two-dimensional swap spectroscopy, a, Optical micrograph of sample, 
showing three coplanar waveguide resonators (R a , R^ and R c , with meander design) capacitively coupled to two super- 
conducting phase qubits (Q 1 and Q 2 , in boxes), b, Block diagram showing main elements, which comprise two circuit 
unit cells R^G^-R^ (light green area) and R^-Q 2 -R c (dark blue area). The horizontal placement represents the spatial 
layout of the sample, whereas the vertical distribution corresponds to the frequencies of the elements. The qubit-resonator 
coupling capacitors are designed to be Cia = ^ib = ^2b = ^2c = 1-9 fF. f a ~ 6.29 GHz, / b ~ 6.82 GHz and f c ~ 6.34 GHz 
are the measured resonator frequencies, and f 1 and f 2 the tuneable qubit transition frequencies. Control and readout 
wiring is also shown, c, Upper sub-panel, pulse sequence for swap spectroscopy, with data shown in d and e. Q 1 line shows 
Gaussian microwave 7r-pulse (red) and qubit tuning pulse (z-pulse; hashed magenta) with variable z-pulse amplitude zpa 
and duration Ar, followed by a triangular measurement pulse (black). Lower sub-panel, diagrammatic representation of 
the pulse sequence. (I) The entire system is initialized in the ground state. (II) Q 1 is excited by a 7r-pulse (red), then (III) 
brought into resonance with R a via a z-pulse (magenta) and allowed to interact with the resonator for a time Ar. (IV) A 
measurement pulse projects the qubit onto its ground state |g) or excited state |e). d, Two-dimensional swap spectroscopy 
for Q 1 . The probability P e to find the qubit in |e) is plotted versus z-pulse amplitude and resonator measurement time 
Ar. The typical chevron pattern generated by a qubit-resonator swap (arrows) is evident for both the Q-^Rg^ (dashed 
black box) and Q x -R^ interactions. Near the center of the plot a qubit interaction with a spurious two-level system is seen, 
surrounded by two regions with short qubit relaxation time, e, Same as d for Q 2 . From these measurements, we find the 
coupling strengths g la ~ 17.58 =F 0.01 MHz, g lh ~ 20.65 =F 0.02 MHz, g 2h ~ 20.43 =F 0.01 MHz and g 2c ~ 17.96 =F 0.01 MHz. 



adjust the qubit operating frequency f 12 and couple microwave pulses for controlling and measuring the 
qubit state. During operation, the device is attached to the mixing chamber of a dilution refrigerator at 
« 25 mK. 

The circuit layout (Fig. |TJo) can be decomposed into two unit cells, Ra-Qi-R^ (g reen area) and R^- 
Q 2 -R c (blue area). The shared resonator connects the two cells and protects the two qubits from 
unwanted crosstalk. The resonator frequencies / a , and f c are measured with qubit spectroscopy (not 
shown). The vacuum Rabi couplings between each qubit and its corresponding resonators, g la and 
for Q x and and g 2c for Q 2 , are determined by their respective coupling capacitors (cf. Fig. |TJd) . The 
coupling strengths are measured using two-dimensional swap spectroscopy (cf. Figjl]i,e). We note that 
the swap spectroscopy provides an excellent tool for revealing the presence of spurious two-level systems 
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Figure 2: Photon shell game, a, Block diagram of the sequence used to coherently transfer a single photon Fock state 
|1) from R a to R c via R^. After initializing the system in the ground state, (I), Q-^ is excited by a 7r-pulse and, (II), 
z-pulsed into resonance with R a for a full Rabi swap (Rabi 7r-swap) at the end of which, (III), R a is populated by the 
one-photon Fock state |1) and Q x is in its ground state, at the idle point. In (IV), Q 1 is z-pulsed into resonance with R a 
for a Rabi 7r-swap and then, (V), z-pulsed into resonance with R^ for another Rabi 7r-swap at the end of which, (VI), R^ 
is populated by the one-photon Fock state and both Q x and Q 2 are in the ground state at the idle point. (VII), Q 2 is 
z-pulsed into resonance with R^ for a Rabi 7r-swap and, (VIII), z-pulsed into resonance with R c for another Rabi 7r-swap 
at the end of which, (IX), R c is in the one-photon Fock state and Q 2 in the ground state at the idle point. Measurement 
and qubit state readout are performed in (X) and (XI), respectively, where the presence of the Fock state in R c is detected 
by its interaction with Q 2 . b, Measurement outcomes for different photon shell games. Each plot shows the probability 
P e to measure a qubit in the excited state |e) as a function of the qubit-resonator measurement time At. Blue circles 
are data, magenta lines a least-squares fit. Data for Q 1 are in the first two columns, for Q 2 in the second two; each row 
corresponds to a different game. Row (i), all resonators are in the vacuum state (with fidelity T > 0.99; we define T as 
the amplitude of the fit; see Supplementary Information), i.e. no stored photons. Row (ii), resonator R a contains one 
photon [T = 0.86 =p 0.01), with the other two resonators in the vacuum state. In row (iii) the photon has been placed in 
R^ (J 7 = 0.80 =p 0.01), and in (iv) the photon is in resonator R c (J 7 = 0.69 =p 0.01). In (v), we have taken the photon from 
R c and placed it back in R^ (J 7 = 0.61 =p 0.01), demonstrating the high degree of control and coherence in the system. 

(TLSs) as well as frequencies with short qubit relaxation times. In all the experiments, the qubits are 
initialized in the ground state |g) and are typically tuned to the idle point, where the qubit Q 1 (Q 2 ) 
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Figure 3: Quantum-mechanical realization of the 'Towers of Hanoi'. Format as in Fig. [2J3, showing the probability 
P e of measuring Q x or Q 2 in the excited state as a function of interaction time At with a resonator. Data are shown as 
blue circles, with a least-squares fit as magenta solid lines, (i), A one-photon Fock state |1) in R a with both and R c 
in the vacuum state, and (ii) a two-photon Fock state in R a (fidelity T = 0.63 =p 0.01) with the other two resonators in 
the vacuum state. A qubit-resonator interaction oscillation is y/2 faster for a two-photon state compared to a one-photon 
state. In (in), one photon has been transferred from R a to R^, so one-photon oscillations are seen when measures R a 
or R^, and when Q 2 measures R^. In (iv), the second photon has been transferred to R^ (J 7 = 0.51 =p 0.01), yielding 
the \[2 increase when either Q-^ or Q 2 measure R^. In (v), both photons have been transferred to R c (T = 0.34 =p 0.01). 
Note that even with this complex protocol, both R a and R^ display negligible oscillations (J 7 > 0.97). See Supplementary 
Information for further analysis. 

|g) O |e) transition frequency f x (/ 2 ) is set in-between, and well away from, the resonator transition 
frequencies / a and f h (/ b and / c ). 

When Q x (Q 2 ) is at the idle point, the qubit-resonator detuning is sufficiently large that the qubit- 
resonator interactions are effectively switched off. A particular qubit-resonator Q p -R q (p = 1,2 and 
q = a, b, c) interaction is switched on by shifting the qubit transition frequency / to equal the resonator 
frequency / , setting the de-tuning to zero and enabling quantum energy transfers. The time-dependent 
control of the qubit transition frequency thus enables highly complex quantum control of the resonators 23 . 

Figure |2^i shows a diagram of the pulse sequence used to implement the single photon equivalent of 
the 'shell game', in which a pea is hidden under one of three shells and the contestant must guess where 
the pea is after the shells have been shuffled. The three resonators play the role of the shells and a single 
photon Fock state |1) that of the pea. The system is initialized in the ground state, with the qubits at 
their idle points, so that all interactions are effectively switched off. Qubit Q x is used to pump a single 
photon into resonator R a [Figure [2^i(I-III)] . The photon state can then be transferred to either of the 
other two resonators, using the qubits in a similar fashion to mediate the single-excitation transfer. A 
transfer from R a to is shown in steps (IV- VI), and a second transfer from R^ to R c performed using 
Q 2 in steps (VII-IX). The final location of the photon can be determined by employing the qubits as 
photon detectors, through the Rabi oscillations that occur when a qubit is brought in resonance with a 
resonator storing a photons [steps (X-XI)]. 

In the data shown in Figure |2b(i)-(v), which represents different versions of the game, a photon was 
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Figure 4: Combined quantum state transfer and storage for one- and two-photon states, a, Probability P e to 
find Q 2 in the excited state (colour bar scale, right side) versus measurement time At (horizontal axis), and total storage 
time T g ^ (vertical axis). A one-photon Fock state is created and stored in resonator R a for a time 7" s ^/3, transferred to 
and stored for the same time, then transferred to R c , stored for the same time and then measured, b, Same as in a, but 
for the generation and storage of a two-photon Fock state. 



stored in one of the three resonators, shuffled between the resonators, and all three resonators then 
measured. In the shell game of Fig. [2j}(i), no photon was placed in any resonator, while for example in 
game (iv), a single photon was transferred from R a to R c via R^, and then detected using qubit Q 2 ; 
measurements of the other resonators R a and R^ show no oscillations, i.e. no photonic excitation. 

We also explored a variant of the shell game, transferring a two-photon Fock state |2) from R a to 
R^ to R c . The two-photon Fock state is first generated in R a 21 , as shown by the measurements in 
Fig. j3p),(ii). Figure |3] also shows the measurements after this state is transferred from R a to R^ and 
then to R c . Each transfer takes two steps, starting with e.g. the state IQiRa^R^) = |g20). The qubit is 
brought into resonance with R a , and one photon is Rabi-swapped to the qubit, at a rate y/2 faster than 
the usual one-photon rate^, leaving the system in the state |elO). The qubit is then tuned into resonance 
with resonator R^ for a one-photon Rabi swap, resulting in the state |gll) [Fig. [3pii)]. The qubit is 
subsequently placed back in resonance with R a for a one-photon swap, yielding |e01), and brought into 
resonance with R^ to transfer the second photon, ending with the state |g02) [Fig. |3|iv)]. To finally 
transfer the photons to resonator R c , the process is repeated using qubit Q 2 , which completes the full 
transfer of Fock state |2) [Fig.^v)]. This process resembles the well-known game 'The towers of Hanoi', 
where a set of disks with different diameters has to be moved between three posts (the three resonators) 
while maintaining the larger disks (Fock state |1), with the longer swapping time) always at the bottom 
of each post, and the smaller disks (Fock state |2), with shorter swapping time) on top. 

A fundamental question for resonator-based quantum computing is whether quantum states can be 
stored in a resonator and later extracted and/or stored elsewhere. We demonstrate this functionality in 
Fig. [4j where a single photon is stored in each of the three resonators for a variable time before 
being transferred to the next resonator. Qubit measurements of the resonator containing the photon 
display clear oscillations for total storage times longer than 3 [is. Figure [4}d is the same experiment 
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Figure 5: Two-resonator Rabi swaps, a, Block diagram of the preparation and measurement protocol. (I), Q 1 is 
excited by a 7r-pulse and, (II), brought into resonance with R a for a Rabi 7r-swap, at the end of which, (III), R a is populated 
by a one-photon Fock state. In (IV), is brought into resonance with R a for a variable transfer time r at the end of 
which, (V), R a is left partially populated, while, (V) and (VI), the remaining energy is fully transferred to R^ via a Rabi 
7r-swap with Q 1 . The probabilities for having the photon in R a or R^ can be varied continuously by changing the transfer 
time r (cf. main text). These probabilities are simultaneously measured for R a with Q ± (VII a) and for R^ with Q 2 (VII 
b). b, Qubit probability P le (colour scale bar) as a function of the Q-^Rg^ measurement time At^ (horizontal axis) and 
the transfer time r (vertical axis), c, Same as b, but for Q 2 's probability P 2e as a function of the Q 2 -R^ measurement 
time Ar 2 . Measurement times are sufficient to display one complete Rabi oscillation between the measurement qubit 
and resonator, with a Rabi swap occurring at the center of each horizontal axis (dashed white line), with multiple swaps 
displayed as a function of the Q 1 -R a transfer time r (vertical direction). Cuts through the probabilities along the dashed 
white lines (full Rabi swaps) are shown in d, where the expected co-sinusoidal oscillations are observed in Q x 's probability 
P le (dark blue circles) (Q 2 's probability P 2e , light green circles), with the summed probability P = P le + P 2e (magenta 
circles) showing the expected slow decay (cf. main text). Solid lines are fits to data. 



repeated with a two-photon Fock state. Clear oscillations are visible for total storage times longer than 
1.5 /is. These experiments demonstrate the realization of a programmable quantum information register. 

We further show that we can generate quantum state entanglement between the two resonators R a 
and R^. The protocol is diagrammed in Fig. |5^l. A one-photon Fock state is first stored in R a , placing 
the system in the state IQiRa^) = |glO) [Fig. [5^(1)- (III)]. The qubit is then used to perform a partial 
transfer of the photon to R^, achieved by placing Q x in resonance with R a and varying the transfer 
time r [Fig. [H^l(IV)]. This leaves the system in the entangled state a|glO) + /3|e00), with amplitudes^ 
a = cos(7r^ la r) and f3 — i sin(7r^ la r). The qubit frequency f x is then tuned from / a to ft , and left there 
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for a time equal to the C^-R^ swap time, thus mapping the qubit state onto the resonator and resulting 
in the two-resonator entangled state 



We then use both qubits to simultaneously measure the two resonators, Q x measuring R a for a 
measurement time Ar t and Q 2 measuring for a time Ar 2 . Figure j5j}(c) shows the resulting oscillations 
in P le (P 2 e) f° r Qi (Q2) (colour bar scale), as a function of the measurement time (horizontal axis) and 
Qi"^a transfer time r (vertical axis). The data display a single Rabi oscillation along the horizontal axis, 
which would repeat if the measurement time were increased, and also shows clear swaps as a function of 
the transfer time r, as expected from the functional dependence of a and f3. If the measurement times 
At-l 2 are chosen to equal a full qubit-resonator swap time (dashed white lines in Fig. J5|d,c), the system 
will be in the state |QiRa^ D Q2) = <^|e00g) + /3|g00e). An ideal measurement of Q 1 and Q 2 would then 
yield probabilities P le = \a\ 2 = cos 2 (7rg la r) and P 2e = \/3\ 2 = sin 2 (7r# la T) for Q 1 and Q 2 , respectively, 
as a function of the transfer time r. In Figure [5]i we display this functional dependence, with a clear 
180° phase difference between the two probabilities and the summed probability P le + P 2e close to unity, 
as expected. The probabilities decrease with r, owing to the finite energy relaxation of Q x and R a . 
The decay time of R^ does not contribute noticeably (cf. Supplementary Information). Fitting yields 
an effective two-resonator decay time approximately equal to the harmonic mean of the relaxation times 



Tf el ~ 340 ns and T| el ~ 3.9 /is of Q 1 and R a , respectively, T 1 ? 1 « (1/27^ + l/2T| e1 )- 1 ~ 626 ns. We 



note that the phase qubit enables a photon transfer between R a and R^ even though the two resonators 
are separated in frequency by ~ 12000 resonator linewidths. 

This last experiment demonstrates a true quantum version of the 'shell game', where the 'pea' (the 
photon Fock state) is simultaneously hidden under two shells, and the contestant's selection of a shell 
constitutes a truly probabilistic measurement. More generally, we have experimentally demonstrated an 
architecture with three resonators and two qubits that displays excellent quantum control over single 
and double microwave photon Fock states. From a fundamental perspective, these results demonstrate 
the potential of multi-resonator circuit QED^ 2 ^, both for scientific study and for quantum information. 



Sample fabrication. The resonators are made from a 150 nm thick rhenium film grown in a molecular-beam-epitaxy 
system onto a polished sapphire substrate. All the other wiring layers are made from sputtered aluminum with A1/A10 X /A1 
Josephson tunnel junctions. All the micro-structures on the different layers are patterned by means of optical lithography 
and etched by means of inductively coupled plasma etching. Amorphous silicon is used as dielectric insulator for the qubit 
shunting capacitors and cross-overs. Our sample fabrication clearly shows the flexibility offered by multi-layer processing. 

The complete device is then bonded by aluminum wire-bonds to an aluminum sample holder, which is bolted to the 

mixing chamber of a dilution refrigerator operating at ~ 25 mK. A detailed description of the fabrication techniques, 

electronics and qubit calibration procedures can be found elsewhereEH. 

Three-resonator circuit QED Hamiltonian. The Hamiltonian H for the circuit of Fig. ^> can be written as the 
sum of the Hamiltonians of each unit cell, H 1 and H 2 - Neglecting the driving and dissipative terms for simplicity, the 
Hamiltonian of the first circuit unit cell can be expressed in the interaction picture with respect to Q 1 and R a and as 
the combination of two Jaynes-Cummings interactions: 



a| g 10> + /3|g01) = |g) ® (a|10> + /?|01» ■ 



METHODS 





(i) 



7 



where &^ are the rising and lowering operators for Q l5 a,a^,b and the bosonic annihilation and creation operators for 
R a and R^, respectively, and A la = f 1 — / a and = f 1 — the qubit-resonator de-tunings. The Hamiltonian H 2 for 
the second circuit unit cell has an analogous expression. 

In order to effectively switch off a particular qubit-resonator Q p -R g interaction, the condition A pq ^> g pq must be 

fulfilled. This is the case when Q 1 (Q 2 ) is at the idle point. On the contrary, when A pq — > 0, a resonant Jaynes-Cummings 

interaction takes place enabling state preparation and transfer in and between the resonators. 
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In this supplementary information, we provide further insight into the transduction fi- 
delity of microwave one- and two-photon Fock states between three coplanar wave guide 
resonators. In addition, we show that the exponential decay of the two-resonator vacuum 
Rabi swaps fits very well with a simple harmonic-mean decay model (confirmed by numer- 
ical simulations), where both the qubit and resonator energy relaxation times contribute 
to the effective decay of the two-resonator dynamics. 

In order to realize three-resonator circuit quantum electrodynamics (QED) experiments, a complex 
architecture must be designed and fabricated on a single chip. In the main text (cf. Fig.[]Ji) we have shown 
an optical micrograph of the main circuit elements, comprising three coplanar wave guide resonators and 
two superconducting phase qubits. Figure [T] shows a detail of qubit Q l5 together with its readout SQUID, 
control and readout lines and the two coupling capacitors to resonators R a and R^. 

In order to unambiguously verify the high fidelity of transfer of a one-photon Fock state from resonator 
R a to resonator R c via resonator R^, we have also performed full-state Wigner tomography on R a and 
R c for the two prototypical examples of photon shell game of Fig. [2jo(ii) ,(iv) (cf. main text). The results 
are displayed in Fig. S.|2j which shows the measured Wigner functions W(a) and corresponding density 
matrices p for the Fock state |1) first stored in R a and then in R c after transfer via R^. 

The Wigner function is obtained following the steps explained in Ref. 1 . In summary, the resonator is 
first prepared in the desired microwave photon state Second, the resonator is displaced by injecting 
a coherent state | — a) with complex amplitude a = \a\ exp(<^ a ), where \a\ represents the coherent state 
real amplitude and cp a its phase; the state is injected through microwave control lines using a classical 
pulsed source (cf. Fig.[l^,b). Third, a qubit in its energy ground state is brought into resonance with the 
resonator for an interaction time long enough to execute several qubit-resonator swaps. Fourth, a simple 
least-squares fit of these oscillations allows us to obtain the probabilities associated with the different 
photon number states, from which it is finally possible to obtain the quasi-probability distributions via 
the parity operator, giving access to full-state Wigner tomography^H The protocols for calibrating the 
amplitude and phase of the coherent state used to displace the resonator are explained in detail in Ref. 
[T| From the Wigner function it is possible to reconstruct the density matrix of the photon state^^ |^). 
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R a Readout 

Fi gure Supplementary 1; Detail of qubit coupled to resonators R a and The capacitors C^ a and 

couple qubit Q x (gradiometer design) to resonators R a and R^. The qubit state is read out by a DC SQUID (gradiometer 
design). A portion of the control and readout lines is also visible. A micrometer scale shows the real circuit dimensions. 

Given the theoretical and measured density matrices of a photon state \Sfr), and p m , respectively, 
we define the state fidelity as T = Tr{p^p m }. For the Fock state |1) prepared in R a , with the 
associated Wigner function and density matrix of Fig. S. [2^i,c, we find T c± 0.84. This compares well 
with the amplitude of the fit to the qubit-resonator swaps, which gives a fidelity T ~ 0.86 (cf. Fig. [| 
in main text). After being transferred to R c , the state has the Wigner function and density matrix of 
Fig. [2]3, d and the corresponding fidelity T ~ 0.64, which is also consistent with the fidelity found with the 
simple least-squares fit (cf. Fig. [2] in main text). The loss of fidelity occurring during the photon transfer 
between the three resonators can be attributed to qubit decoherence and slight calibration errors during 
the swap qubit-resonator operations. Nevertheless, it is remarkable that the density matrix associated 
with the state in R c is still very pure, with nearly negligible spurious matrix elements and only a small 
contribution from the |0) state. 

Other two relevant figures of merit for the transfer of one- and two-photon Fock states between 
resonators R a , R^ and R c are represented by the harmonic purity of the state (i.e. absence of beatings) 
and the y/2 scaling factor between swaps associated with the Fock states^ |1) and |2). These properties 
can easily be estimated by computing the fast Fourier transform (FFT) of the qubit-resonator vacuum 
Rabi swaps. Figure Supplementary [3] shows the FFTs of the Rabi swaps for the key steps of the photon 
shell game and 'Towers of Hanoi', i.e. for a Fock |1) and |2) created in R a and measured with Q l5 then 
transferred to R^ and measured with Q x and Q 2 , and finally transferred to R c and measured with Q 2 . 
The time-domain swaps used to compute the FFTs are those in Figs.|2]3 and [3] in main text. 
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Figure Supplementary 2: Wigner tomography for the photon shell game, a, Measured Wigner function W(a) 
for resonator R a as a function of the complex resonator amplitude a in square root of photon number units (colour scale 
bar on the far right). The inset displays a cut of the three-dimensional plot of the Wigner function, b, Same as in a, 
but for resonator R c (colour scale bar on the right). Negative quasi-probabilities clearly indicate the quantum-mechanical 
nature of the intra-resonator states, c, Theoretical (grey) and measured (black) values for the density matrix associated 
with the state stored in resonator R a , p, projected onto the number states p mn = (m\p\n). The magnitude and phase of 
p mn is represented by the length and direction of an arrow in the complex plane (the scale for the real and imaginary part 
is reported on the far right), d, Same as in c, but for resonator R c . When representing the density matrices, the resonator 
Hilbert space has been truncated to the lowest four bosonic states. 



The last topic to be considered in this supplementary is the energy relaxation mechanism of the 
two-resonator Rabi swaps of Fig. [5]i in the main text. 

In order to gain deeper insight into the decay dynamics of the two-resonator Rabi, we first need infor- 
mation on the energy relaxation times of qubit Q x and resonators R a and R^. Figure Supplementary |4j 
shows the experimental measurement of these decay times. Note that the resonators' energy relaxation 
time is determined by preparing the resonator in a one-photon Fock state |1), and storing it for a variable 
time. After this time, a qubit is brought on resonance with the resonator, the state swapped into the 
qubit, and the qubit population then read out^l 

The energy relaxation times obtained from a simple exponential fit are reported in Table S. [IJ We 
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Q r R a Q r R b Q 2 -R b Q 2 -R c 




Vacuum Rabi frequency (MHz) 

Figure Supplementary 3: Fourier analysis for the photon shell game and 'Towers of Hanoi'. Top sub-panels, 
normalized Fourier amplitude as a function of vacuum Rabi frequency associated with the qubit-resonator Rabi swaps for 
a one-photon Fock state |1). Bottom sub-panels, normalized Fourier amplitude for a two-photon Fock state |2). Above 
each column is indicated the respective qubit-resonator interaction. The dashed black line in each sub-panel indicates the 
maximum Fourier component. The amplitude re-normalization is calculated with respect to the largest Fourier component 
for Fock state |1) (top sub-panels) and \2) (bottom sub-panels), respectively. The \f2 scaling^ between the |1) and |2) 
states is clearly visible. The absence of beatings (only a small beating, owing to some residual presence of state |1), for state 
1 2) in the Q 2 -R^ sub-panel) shows the high level of harmonic purity of the states transferred between the three resonators, 
both for the photon shell game and for the more complex 'Towers of Hanoi'. 



use the qubit and resonator energy relaxation times as well as the other parameters listed in the table 
to solve numerically a Lindblad-type master equation^^l: 

1 ^ ^ 3 

P=^( H iP-P H i)^J2^P, (2) 

k=l 

where p = (d/dt) p is the time derivative of the total density matrix p describing the system, H 1 is the 
Hamiltonian of Eq. ([!]) (cf. Methods' section in main text), t k is the Lindblad superoperator defined 
as C k p = j k (X k pXl — X k X k p/2 — pX\X k /2) and k G N. The qubit and resonators decay rates are 
defined as 7 X = 1/T| e ^, 7 2 = l/T-p^ and 7 3 = and the generating operators as X 1 = a, X\ = a\ 

X 2 = a~ , X.\ = <t+ , X 3 = b and X\ = P. We numerically solve Eq. ^ for the pulse sequence shown 
in Fig. [5^i (cf. main text), without accounting for the measurement process. The results are displayed in 
Fig. S.[5| where they are compared to the experimental data. 

Figure Supplementary [5^i,c shows the same data as in Fig. |5]i, but for a transfer time r three times 
longer. Furthermore, the data shown in Fig. S. [5^i,c were taken using a different device compared to 
Fig. [5]i, with longer qubit relaxation times. The exponential decay obtained by the simple harmonic 
mean model (cf. main text) is superposed with the data, making evident the qualitative validity of the 
model. Figure Supplementary [5]a,d shows the results of the simulations of Eq. Q corresponding to 
the experimental data of Fig. S. |5^i,c, respectively. Data and simulations are in very good agreement, 
supporting the simple harmonic mean decay model. In particular, the experimental decay time obtained 
by fitting the data is ~ 840 ns, from simulations ~ 874 ns and from the harmonic mean model ~ 896 ns. 
The amplitude of the simulated two-resonator Rabi swaps was adjusted to the measured amplitudes. The 
slight discrepancy between the experimental data and simulations for the low region of the occupation 
probabilities (causing an offset between data and simulations) is because the simulations do not account 
for the measurement process. Notice that we can safely assume that only qubit Q x and resonator R a , 
swapping for a variable transfer time r, contribute to the effective decay mechanism of the two-resonator 



14 




Fi gure Supplementary 4: Energy relaxation for qubit Q 1 and resonators R a and In all panels: P pe is 

the probability to find Q p in |e) as a function of measurement delay time (i.e. the time a photon is stored in the resonator 
before being swapped into the qubit). Full circles are data and solid magenta lines exponential fits to data, a, Measurement 
of the energy relaxation of resonator R a using qubit Q x as a detectoi>^. b, Energy relaxation of qubit Q 1 . c, Same as in 
a, but for resonator R, . 



Rabi dynamics. In fact, the second resonator serves only as a mapping resonator, the state of which 
is measured with Q 2 typically after a time Ar 2 <C (cf. main text and Fig. j^c). In other words, 
examining the two-dimensional plots of Fig. [5}d,c we expect two distinct decay mechanisms. The first 
along the horizontal axis (Ar x and Ar 2 ). This decay is practically negligible as this measurement is 
completed in ^ 30 ns. The second is along the vertical axis (r) related to the Qi-Ra swaps, as explained 
above. 

For a theoretical analysis of the decay mechanisms characteristic of two-resonator dynamics in differ- 
ent regimes we remind to Ref. [I3j 
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Table Supplementary 1: Main parameters used for the numerical simulations of the two-resonator Rabi 
dynamics. S represents the non-linearity of the superconducting phase qubit, i.e. the energy difference in unit Hz between 
the qubit transition frequencies relative to the ground state and first excited state and to the first excited state and second 
excited state, respectively^. The non-linearity has been used in the simulations to take into account possible leakage 
outside of the qubit subspace. and are the dephasing times for resonators R a and R^, respectively. Since our 

main goal is to understand the energy relaxation of the two-resonator Rabi swaps, the qubit dephasing time Tf has been 
neglected in the simulations. All the other parameters are defined in the main text. 
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As a last remark, it is worth mentioning that all data shown in the main text and supplementary 
information were corrected for measurement errors, following the standard procedure outlined in Ref.fTTl 

In conclusion, in this supplementary information we have validated the experimental results and 
conclusions presented in the main text by means of full-state Wigner tomography, FFT analysis and 
numerical simulations. 
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Fi gure Supplementary 5; Comparison between experimental results and numerical simulations of two- 
resonator Rabi swaps, a, R a state dynamics measured as the probability P le to find Q x in the excited state versus 
transfer time r. Dark blue circles are data, solid magenta line a least-squares fit to data and dashed black line the 
exponential harmonic mean decay, b, Simulation (parameters in Table S.Q of the data in a, showing the resonator mean 
photon number (a^a) versus r. The exponential decay envelope (dashed black line) confirms the harmonic mean decay 
model, c, Same as in a, but for (cf. main text). Light gr^en circles are data, solid magenta line a least-squares fit to 
data and dashed black line the exponential harmonic mean aecay. d, Same as in b, but for R^. 
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